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ABSTRACT 

A  primitive-equation.  n-layer,  eddy-resolving  circulation  model  has  been  applied  to  the  (lull  Stream  System 
from  (ape  llatleras  to  east  of  the  ( irand  Banks  ( 78°-45"W.  10” -4X"  N  ).  Within  the  limitations  of  the  model, 
realistic  coastlines,  bottom  topography ,  and  forcing  functions  have  been  used.  A  two-layer  version  of  the  model 
was  driven  by  observed  mean  climatological  wind  forcing  and  mass  transport  prescribed  at  inflow.  Outflow  was 
determined  by  a  radiation  boundary  condition  and  an  integral  constiaint  on  the  mass  field  in  each  layer 
Specification  of  a  IVep  Western  Boundary  Current  (I)WBC)  was  included  in  some  model  runs. 

Six  numerical  experiments,  from  a  series  of  over  fifty  integrated  to  statistical  equilibrium,  were  selected  for 
detailed  description  and  intercomparison  with  observations.  The  basie  ease  consisted  of  a  flat  bottom  regime 
driven  by  wind  forcing  only  Realistic  inflow  transport  in  the  upper  layer  was  then  prescribed  and  two  different 
outflow  specifications  at  the  eastern  boundary  wen-  studied  in  experiments  2  and  .V  Three  additional  experiments 
involved  < 4 1  adding  bottom  topography  ( ineludi  v  the  New  Hngland  Seamount  (  hain ),  (  5 )  adding  a  DWB( ' 
to  experiment  4  with  20  Sv  (Sv  »  10'’  m  ’  s  1 )  total  transport,  and  (6)  increasing  the  DWBC  to  40  Sv.  A  brief 
discussion  of  the  influence  of  parameter  variations  includes  modifications  of  dissipation  (lateral  eddy  diffusion 
and  bottom  friction)  and  stratification. 

Results  fiom  the  sequence  of  experiments  suggest  an  important  role  for  the  DWBC  in  determining  the  mean 
path  of  the  <  lull  Stream  and  consequently  the  distribution  of  eddy  kinetic  energy,  and  the  character  of  the  deep 
mean  flow.  The  most  realistic  experiment  compares  to  within  a  factor  of  two  or  better  with  observations  of  the 
amplitude  of  eddy  kinetic  energy  and  rms  fluctuations  of  the  thermocline  and  sea  surface  height.  Abyssal  eddy 
kinetic  energ-  was  smaller  than  observed.  The  mean  flow  is  characterized  by  recirculations  to  the  north  and 
south  of  the  Gulf  Stream  and  a  deep  cyclonic  gyre  just  cast  of  the  norther.  portion  of  the  New  Hngland 
Seamount  Chain,  as  found  in  the  data. 


I.  Introduction 

The  dynamics  associated  with  instabilities  of  western 
boundary  currents  and  the  role  of  the  resultant  mc- 
soscale  eddies  in  understanding  the  ocean  general  cir¬ 
culation  are  important  topics  in  modern  oceanography. 
Recently,  eddy-resolving  models  have  become  suffi¬ 
ciently  realistic  and  the  data  sufficiently  extensive  so 
that  modellers  and  observationalists  have  begun  to 
compare  their  respective  results,  notably  those  relevant 
to  the  Gulf  Stream  System.  In  a  series  of  studies  (e.g., 
Schmitz  and  Holland  1982;  Holland  and  Schmitz.  1985, 
Schmitz  and  Holland  1986;  henceforth  designated 
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SH82,  HS85,  and  SH86,  respectively)  results  from 
mostly  rectangular-domain  quasigeostrophic  numerical 
models  were  compared  with  observations  of  time  av¬ 
erages  from  the  western  North  Atlantic  and  Pacific.  In 
essence  these  papers  addressed  the  question:  Can  an 
eddy-resolving  wean  model  simulate  the  observed  dis¬ 
tributions  of  ocean  variability  as  measured,  for  ex¬ 
ample,  by  eddy  kinetic  energy  amplitudes,  mean  flows, 
and  distributions  of  the  off-diagonal  components  of 
the  horizontal  Reynolds’  stresses? 

The  eddy  kinetic  energy  distributions  in  the  mid¬ 
latitude  North  Atlantic  are  now  known  in  general  out¬ 
line.  for  example,  fig.  la  shows  eddy  kinetic  energy 
( f’Kf )  calculated  from  surface  drifters  (Richardson 
1983).  At  the  surface  there  is  a  small  region  of  maxi¬ 
mum  kinetic  energy  greater  than  2000  cm:  s  :  (and 
perhaps  greater  than  3000  cm :  s  : )  west  of  the  New 
F.ngland  Seamount  Chain  (NKSC).  These  results  arc 
in  rough  agreement  with  estimates  of  F.Kf.  from  alti- 
metric  measurements  (GFOS-3)  by  Douglas  et  al. 
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Fir..  1.  (a)  Eddy  kinetic  energy  of  the  upper  ocean  (in  cm-’  s  ’)  determined  from  drifters  (  Richardson  14SJ).  (hi 
Updated  chart  of  the  existing  estimates  of  abyssal  eddy  kinetic  energy  (cm;  s  i*  the  North  Atlantic.  (See  also 
Schmitz  1984). 


( 1983).  but  are  higher  than  those  estimated  from  ship 
drift  ( Wyrtki  et  al.  1976).  Figure  1  b  is  an  updated  map 
of  abyssal  eddy  kinetic  energy  estimates  for  the  North 
Atlantic.  Most  of  the  data  sources  have  been  described 
in  Schmitz  (  1984).  The  maximum  of  140  ctrr  s  :  in 
EKE  east  of  the  NESC  at  55°W  at  the  latitude  of  the 
mean  Gulf  Stream  is  notable,  as  is  a  two  order  of  mag¬ 
nitude  latitudinal  decay  into  the  interior  of  the  sub¬ 


tropical  gyre  in  both  zonal  and  meridional  directions. 
In  addition  to  the  information  on  EKE.  a  coherent 
picture  of  the  deep  mean  flow  in  the  western  North 
Atlantic  is  beginning  to  emerge.  For  example.  Fig.  Ic 
is  a  map  of  long-term  current  meter  measurements 
from  Hogg  et  al.  ( 1986).  Also  shown  is  the  mean  po¬ 
sition  of  the  15°C  isotherm  at  200  m  adapted  from 
Fisher  (  1977).  The  DWBC  flowing  along  the  conti- 
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cm/sec 

IT< I .  ( (  ontnun't/)  u')  Deep  mean  velocities  from  measurements  compiled  by  Hogg  et  al.  (  19X6 ).  The  linn  portion 
of  the  arrow  is  an  estimate  ol  the  uncertainly  in  the  mean  resulting  from  eddy  noise.  Also  shown  (solid  thin  line)  is 
the  mean  position  of  the  1 5“C.  isotheim  at  200  m  adapted  from  l  isher  ( 1977  ).  The  isotherm  is  between  the  dashed 
lines  50%  of  the  time 


ncntal  rise,  an  offshore  “northern  recirculating  gyre.”  were  different  in  shape,  possibly  due  to  the  New 
and  a  strong  westward  How  north  of  about  40°N  just  England  Seamount  Chain  or  topography  in  general, 
east  of  the  Seamount  Chain  with  no  counterpart  at  Model  equivalent  isotherm  depth  fluctuations  were 
similar  depths  further  to  the  west  are  all  noted  by  Hogg  smaller  than  observed. 

el  al.  (  1986).  An  anticyclonic  circulation  around  the  The  selection  of  precise  geographical  locations  at 
northeast  extension  of  the  Bermuda  Rise  centered  near  which  to  compare  model  and  observations  has  not  been 
35°N,  60° W  was  also  identified  in  the  data.  straightforward  since  model  geometry  and  wind  forcing 

In  general  the  early  two-layer  model  solutions  of  have  been  idealized.  Until  recently  most  of  the  mod- 
SII82  possessed  the  basic  property  that  the  most  encr-  cling  experiments  have  utilized  quasigeostrophic  dv- 
getic  eddies  were  located  near  strong  curicnts  which,  namics  in  regular  domains  with  idealized  winds,  with 
in  turn,  were  fed  by  return  circulations.  However,  the  the  mean  path  of  the  modeled  Gulf  Stream  being  en- 
zonal  scale  for  both  the  mean  and  the  eddy  How  were  tirelv  zonal,  coincident  with  the  zero  curl  line  of  the 
smaller  than  observed.  Hie  HS8?  experiments  dom-  mean  wind  forcing.  Comparison  with  actual  wind  fields 
onstrated  that  increased  vertical  resolution,  increased  in  the  North  Atlantic  suggests  that,  even  with  the  un- 
lateral  friction,  and  introduction  of  a  gentle  bottom  certainties  in  the  data,  the  Gulf  Stream  mean  path  is 
slope  can  improve  the  penetration  scale  of  the  model  significantly  different  from  the  zero  curl  line  of  the 
Gulf  Stream  by  confining  the  input  of  wind  energy  to  mean  annual  wind  and  distinctively  nonzonal.  Also,  a 
shallower  upper  layers  and  inhibiting  instability  pro-  condition  of  the  QG  approximation  is  that  bottom  to- 
cesses.  In  the  higher  resolution  (eight  layer)  experi-  pography  must  be  of  small  amplitude,  no  more  than 
ments  of  SIIK6.  zonal  and  vertical  distributions  of  eddy  10";  of  the  fluid  depth.  Such  a  condition  is  inappro- 
kinclic  energy  near  both  the  Kuroshio  Extension  and  priate  near  the  continental  slope  or  within  (he  New 
the  Gulf  Stream  were  in  reasonable  agreement  with  England  Seamount  Chain  ( NESC),  where  seamounts 
observations,  though  better  for  the  North  Pacific  than  rise  3  km  above  the  surrounding  sea  (loor  and  come 
the  North  Atlantic.  However,  mean  zonal  currents  at  within  I  km  of  the  sea  surface,  furthermore.  Harrison 
some  locations  in  the  midlatitude  jet  were  clearly  higher  (  1982)  has  shown  that  there  can  be  no  net  vorticity 
than  observed  and  the  data-based  surface  KE  distri-  input  to  the  lower  layer  of  a  QG  model  via  vortex 
bution  along  the  Gulf  Stream  and  the  model  analogue  stretching. 
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Thompson  and  Hurlburt  ( 1982)  and  Hurlburt  and 
Thompson  ( 1984),  henceforth  referred  to  as  TH82  and 
HT84,  developed  a  primitive  equation  model  of  the 
Gulf  Stream  including  large-amplitude  bottom  topog¬ 
raphy.  In  the  HT84  study  a  rectangular  model  domain 
rotated  28  degrees  counterclockwise  from  zonal  ex¬ 
tended  from  Cape  Hatteras  to  the  Grand  Banks.  The 
model  domain  was  closed  except  for  an  inflow  port  at 
Cape  Hatteras  and  an  outflow  port  roughly  at  the  mean 
position  of  the  Gulf  Stream  as  it  passes  south  of  the 
Grand  Banks,  in  all  of  the  experiments  the  flow  of  the 
Gulf  Stream  was  confined  to  the  upper  layer  of  the 
model  and  did  not  directly  impinge  on  the  topography. 
Flow  in  the  lower  layer  was  driven  from  above  by  pres¬ 
sure  fluctuations  associated  with  the  active  upper  ocean. 
The  model  NESC  had  a  substantial  influence  on  the 
Gulf  Stream  variability  and  mean  path.  Results  from 
the  reduced-gravity  (equivalent  barotropic)  model  were 
unrealistic,  with  length  scales  for  rings  and  meanders 
too  large  by  a  factor  of  two.  Addition  of  the  barotropic 
mode  and  the  possibility  for  baroclinic  instability 
yielded  a  more  realistic  GS.  Warm  core  ring  generation 
was  enhanced  by  the  NESC,  with  prolific  generation 
occurring  immediately  upstream  of  the  NESC  with  a 
frequency  of  about  2  yr  '.  When  compared  to  the  flat 
bottom  experiment,  the  NESC  increased  both  baro¬ 
clinic  and  barotropic  instabilities  as  measured  by  the 
rate  of  transfer  of  mean  potential  to  eddy  potential 
energy  and  mean  kinetic  to  eddy  kinetic  energy. 

Maximum  EKE  in  the  large  seamount  experiment 
exceeded  2000  cm 2  s  ~2  in  the  upper  layer  and  1 30  cm 2 
s“2  in  the  lower  layer.  Maximum  sea  height  variability 
greater  than  35  cm  was  found  in  the  large  amplitude 
seamount  experiment,  comparable  to  that  observed  by 
Douglas  et  al.  ( 1983).  These  values  are  comparable  to 
observations  but  location  of  the  abyssal  maximum  was 
upstream  of  the  seamounts,  unlike  the  observations, 
and  the  flow  east  of  the  NESC  was  generally  much  too 
weak.  An  unresolved  question  in  HT84  was  the  influ¬ 
ence  of  the  outflow  boundary  condition  on  the  Gulf 
Stream  mean  path  and  energetics.  Because  the  outflow 
was  confined  to  a  narrow  port  at  approximately  the 
latitude  of  the  mean  Gulf  Stream  it  was  possible  that 
the  mean  path  and  the  meanders  of  the  Stream  (and 
thus  KE  east  of  the  NESC )  were  artificially  constrained. 

Adamec  ( 1988)  examined  the  effects  of  seamounts 
in  a  two-mode  quasigeostrophic  model  and  found  that 
seamounts  deflected  the  average  model  jet  southward 
and  induced  a  standing  wave  downstream.  Additional 
vertical  modes  enhanced  the  baroclinic  instability 
mechanism  and  reduced  the  eastward  penetration  of 
the  jet.  However,  in  contrast  to  HT84,  the  seamounts 
in  the  quasigeostrophic  model  had  little  effect  on  the 
distribution  and  magnitude  of  the  eddy  kinetic  energy 
at  the  surface. 

This  paper  is  a  logical  extension  of  HT84,  SH82. 
and  HS85.  Our  purpose  is  to  describe  the  experimental 


design  and  some  results  from  a  new  primitive-equation 
regional  model  of  the  Gulf  Stream  including  realistic 
geometry,  bottom  topography,  and  wind  forcing.  The 
model  results  are  unique  in  that  they  represent  statis¬ 
tically-equilibrated,  eddy-resolving  integrations  of  an 
open  ocean  region  in  an  irregular  domain.  Although 
analysis  and  experimentation  are  ongoing,  we  report 
our  initial  efforts  in  this  paper  as  follows.  Section  2 
describes  the  model  domain,  forcing,  parameters,  and 
inflow  and  outflow  boundary  specifications.  Section  3 
provides  results  from  six  experiments  from  a  series  of 
over  fifty  and  compares  model  results  with  observa¬ 
tions.  focusing  on  areas  of  significant  agreement  and 
disagreement.  One  basic  experiment  serves  as  a  bench¬ 
mark  about  which  bottom  topography,  a  deep  western 
boundary'  current,  and  modifications  to  the  extent  of 
the  open  outflow  boundary  condition  are  added.  In 
section  4  we  summarize  our  results,  identify  deficien¬ 
cies  in  our  approach,  and  suggest  some  future  work. 
Finally,  appendix  B  describes  the  numerical  design  and 
the  open  boundary  conditions  in  some  detail. 

2.  Experimental  design 

The  HT84  Gulf  Stream  model  was  adapted  from  an 
earlier  rectangular-domain  primitive-equation  for¬ 
mulation  for  the  Gulf  of  Mexico  by  Hurlburt  and 
Thompson  (1980.  hereafter  HT80).  Wallcraft  and 
Thompson  ( 1984)  extended  the  model  to  include  ir¬ 
regular  coastlines,  yielding  the  first  eddy-resolving 
primitive-equation  ocean  model  with  realistic  geom¬ 
etry.  bottom  topography  and  open  boundary  condi¬ 
tions  capable  of  multiyear  integrations  to  statistical 
equilibrium.  Benchmark  experiments  using  that  model 
have  been  described  by  Thompson  (1986)  and  by 
Hurlburt  ( 1984). 

Differences  between  the  present  model  and  Hurlburt 
and  Thompson  ( 1980)  are  described  in  appendix  R 
The  model  domain  is  shown  in  Fig.  2.  Inflow  and  out¬ 
flow  specifications  are  described  in  section  3.  Since  a 
focus  of  this  study  is  the  influence  of  the  DWBC  on 
the  Gulf  Stream,  we  have  modeled  it  via  a  specified 
transport  in  the  deep  layer  centered  along  the  3000  m 
isobath  of  the  continental  slope  in  the  northeast  portion 
of  the  model  with  outflow  along  the  southern  boundary . 
Although  only  two-layer  results  will  be  reported  here, 
the  model  has  been  formulated  for  arbitrary  vertical 
resolution  and  three-layer  experiments  are  underway  . 
As  is  usual  with  adiabatic  models  with  a  Lagrangian 
vertical  coordinate,  tonography  is  confined  to  the 
deepest  layer  of  the  system,  and  ventilation  of  the 
model  thcrmocline  is  not  permitted.  If  the  upper  layer 
thins  to  zero  thickness  the  experiment  must  be  ter¬ 
minated.  It  is  possible  to  prevent  such  behavior  and 
continue  integrations,  even  in  an  adiabatic  model,  by 
employing  techniques  in  shock  wave  modeling  such  as 
the  Flux  Corrected  Transport  Method  ( Boris  and  Book 
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Fig.  2.  Model  domain  for  the  limited-area  Gulf  Stream  model.  Bottom  depths  are  shown  with 
a  200  meter  contour  interval.  Inflow  ( T,in)  and  outflow  ( T,ou' )  location  for  upper  ( i  =  1 )  and 
lower  (i  =  2)  layer  transports  are  also  noted. 


1973).  The  method  has  been  applied  by  Bogue  et  al. 
( 1986)  for  the  problem  of  thermocline  outcropping  of 
the  Gulf  Stream  but  we  have  sidestepped  the  problem 
in  these  experiments  by  choosing  a  sufficiently  thick 
upper  layer  at  the  initial  time  to  avoid  thermocline 
surfacing. 

For  the  experiments  with  bottom  topography  the 
original  SYNBAPS  bathymetric  data  base  (Van 
Wyckhouse  1979)  was  mapped  to  a  0.2  degree  model 
grid.  The  model  lateral  boundary  was  taken  to  be  lo¬ 
cated  at  the  200  m  contour.  Except  at  the  ports  the 
boundaries  are  rigid  and  the  no-slip  condition  was  im¬ 
posed.  A  five-point  smoother  was  applied  twice  to  the 
full  topography  to  damp  high  spatial  variability  and 
the  result  scaled  such  that  the  shallowest  topography 
in  the  basin  is  at  1500  m  (see  Fig.  2). 

Model  parameters  are  provided  in  Table  !  and  a 
summary  of  the  numerical  experiments  are  listed  in 
Table  2.  While  it  is  not  our  purpose  to  conduct  an 
exploration  of  parameter  space  in  this  paper,  it  is  im¬ 
portant  to  mention  some  of  the  model’s  sensitivities 
to  parameter  variations.  First,  the  value  of  the  density 
contrast  between  layers  is  based  on  the  vertical  density 
variation  for  a  rather  arbitrary  mean  density  profile 
near  the  “center  of  the  recirculation  gyre.”  If  the  density 


Table  1.  Model  parameters. 


A 

10f  cm2  s  1 

K 

980  cm  s  2 

S' 

2  cm  s  2 

//, 

1000  m 

Hi 

4000  m 

Xx.  Xy 

0.2° 

Xt 

(i0  min 

C„ 

0  and  2x10' 

contrast  is  reduced  the  thermocline  will  eventually  sur¬ 
face.  If  the  density  contrast  is  increased  the  thermocline 
fluctuations  become  too  smali  and  unrealistic.  With 
respect  to  horizontal  eddy  viscosity  for  the  Laplacian 
lateral  friction,  we  have  chosen  the  lowest  value  ( 100 
m  s“2)  which  was  numerically  stable  for  the  0.2  degree 
horizontal  resolution.  The  maximum  grid  interval 
Reynolds  number  for  this  choice  of  lateral  friction  is 
around  1 50.  For  the  quadratic  bottom  friction  we  have 
used  a  stress  coefficient  (C*)  ranging  from  zero  to  2 
X  10“3.  In  the  flat  bottom  experiments  it  is  necessary 
to  include  bottom  friction  to  insure  statistical  equilib¬ 
rium  at  realistic  energy  levels.  Otherwise  baroclinic  in¬ 
stabilities  transfer  energy  to  the  deep  water  at  a  rate 
which  cannot  be  dissipated  via  lateral  friction  or  ex¬ 
ported  from  the  domain.  However,  in  the  experiments 
with  bottom  topography,  baroclinic  instability  pro¬ 
cesses  tend  to  be  inhibited  and  energy  transfer  to  the 
deep  water  is  slowed.  Consequently  bottom  friction  is 
not  necessary  for  statistical  equilibrium  in  these  ex¬ 
periments.  Experiments  5a  (C*  =  0)  and  5b  (Q  =  2 
X  103)  illustrate  this  model  sensitivity. 

3.  Experimental  results 

The  extent  to  which  local  wind  forcing  alone  deter¬ 
mines  the  mean  path,  strength  of  the  Stream  and  its 
recirculation,  and  the  intensity  of  the  eddy  activity  in 
the  absence  of  forced  inflow  is  the  focus  of  Experiment 
1 .  The  model  was  driven  from  rest  with  mean  annual 
wind  forcing  derived  from  12-hourly  FNOC  surface 
(19m)  analyses  from  1 979  to  1 984.  The  resulting  wind 
stress  curl  field  is  shown  in  Fig.  3a.  The  zero-curl  con¬ 
tour  is  nonzonal  and  the  wind  curl  distribution  is  sim¬ 
ilar  to  that  derived  by  Hellerman  and  Rosenstein 
(1983). 
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Table  2.  Model  experiments. 

Volume  19 

Experiment 

c„ 

Transport  (Sv)* 

(Layer  1  /Layer  2) 

Topography 

Port  forcing 

1 

2  X  10'3 

0/0;  wind  forced 

No 

closed  boundaries 

2 

2  x  i<r3 

50/0 

No 

narrow  outflow 

3 

2  X  to-3 

50/0 

No 

wide  outflow 

4 

0 

50/0 

Yes 

wide  outflow 

5a 

0 

50/20 

Yes 

wide  outflow 

5b 

2  X  10-3 

50/20 

Yes 

wide  outflow 

6 

0 

50/40 

Yes 

wide  outflow 

Transport  profiles  on  inflow  are  parabolic. 

Lower  layer  transport  represents  the  Deep  Western  Boundary  Current  flowing  along  the  continental  slope. 
*  Transport  on  Inflow  =  Transport  on  Outflow  for  each  layer. 


Figure  3b  shows  the  synoptic  sea  surface  height  served  (Schmitz  1984),  having  maximum  values  less 
(SSH)  for  Experiment  1  following  spinup  to  statistical  than  11  cm2  s-2  (compared  to  140  cm2  s~2  from  ob- 
equilibrium  (requiring  about  6  years).  To  first  order  servations),  and  is  confined  to  west  of  60°W. 
this  pattern  is  also  the  pattern  of  model  thermocline  Experiment  2  utilizes  the  addition  of  the  inflow 
anomaly,  with  a  10  cm  change  in  SSH  equivalent  to  transport  south  of  Cape  Hatteras.  The  details  of  the 
about  a  50  m  change  in  thermocline  anomaly.  There  inflow  and  outflow  boundary  conditions  are  in  appen- 
isa  weak  Gulf  Stream  [about  20  Sv  (Sv  =  106  m3  s'1 )  dix  B.  The  model  supports  specification  of  either  trans- 
on  average,  consistent  with  the  mean  wind  curl]  ex-  port  or  current  profile  plus  the  inflow  angle  and  requires 
tending  from  Cape  Hatteras  northeastward  and  intense  conservation  of  mass  in  the  domain  at  all  times.  We 
recirculations  north  of  the  stream  near  the  coast  and  specify  a  steady  inflow  transport  and  cross-stream  pro¬ 
in  the  southwest  corner  of  the  domain.  Instantaneous  file  based  on  data  from  Richardson  et  al.  ( 1969).  Total 
penetration  of  the  Stream  into  the  domain  only  extends  upper  layer  inflow  transport  is  50  Sv,  distributed  over 
to  about  65°W.  Eddy  activity  is  evident  from  the  SSH  a  100  km  wide  inflow  with  a  parabolic  cross-stream 
fieid  in  the  weak  meandering  of  the  Stream  and  in  the  profile,  and  an  inflow  angle  of  28  degrees.  The  model 
cyclonic  eddy  centered  al  about  34°N,  73° W.  The  eddy  is  spun  up  with  an  e-folding  time  of  one  year  (see  ap- 
field  is  much  more  clearly  shown  in  the  subthermocline  pendix  B).  Surprisingly,  the  model  was  not  very  sen- 
pressure  field  (Fig.  3c).  The  basin  is  nearly  filled  with  sitive  to  inflow  angle  (although  more  experiments  are 
eddies.  Although  a  large  barotropic  Rossby  wave  re-  necessary  to  understand  inflow  angle  dependencies  and 
sponse  is  excited  during  the  early  spinup,  the  deep  eddy  sensitivities  to  time  dependent  inflow  transports ) . 
field  shown  here  in  year  9  is  energized  via  instabilities  Two  variations  of  the  outflow  condition  were  con- 
of  the  Gulf  Stream  and  associated  wave  radiation.  sidered.  The  first,  used  in  Experiment  2.  prescribed  a 
The  mean  SSH  is  shown  in  Fig.  4a.  The  tight  recir-  narrow  ( 100  km  wide)  outflow  boundary  port  centered 
culation  to  the  south  of  Cape  Hatteras  is  due,  in  part,  at  39°N,  approximately  the  location  of  the  elimato- 
to  the  influence  of  the  rigid  wall  at  30°N.  An  important  logical  mean  Gulf  Stream  at  45  °W.  This  is  similar  to 
aspect  of  the  mean  flow  is  the  cyclonic  recirculation  the  approach  used  earlier  in  HT84.  The  second  con- 
to  the  north  of  the  Stream  and  the  southward  flow  dition  allows  the  entire  eastern  boundary  in  the  upper 
along  the  coast  between  Cape  Hatteras  and  Cape  Cod.  layer  to  be  open  and  is  used  in  the  remainder  of  the 
As  will  be  shown  later,  a  critical  determinant  for  sep-  experiments.  The  requirements  for  mass  balance  and 
aration  of  the  Gulf  Stream  is  the  intensity  of  the  south-  the  possibility  of  a  return  circulation  are  discussed  in 
wr_1  flow  in  this  region,  doth  at  the  surface  and  below  appendix  B. 

the  thermocline  (in  the  DWBC).  Experiment  2  was  integrated  to  statistical  equilib- 

Maps  of  the  eddy  kinetic  energy  (  EKE )  above  and  rium  as  before,  f  igure  5a  shows  a  snapshot  of  the  sea 
below  the  thermocline  (upper  and  lower  layers)  for  surface  height  in  model  year  20.  The  meanders  of  the 
Experiment  I  are  shown  in  Fig.  4b  and  4c.  The  upper  Stream,  warm  and  cold  rings,  and  the  recirculation  are 
layer  EKE  is  very  weak,  about  10G  of  the  observed  evident.  The  maximum  SSH  variation  is  over  140  cm 
surlace  EKE  ( Richardson  1983 )  and  fails  to  penetrate  and  the  height  change  across  the  Stream  is  comparable 
east  of  60° W.  (Note  the  model  “surface"  values  are  to  observation,  about  one  meter.  The  mean  SSH  in 
really  a  layer  average  over  1000  m  in  the  mean  and  Fig.  5b  shows  a  standing  wave  for  the  Stream  with  a 
thus  cannot  be  directly  compared  with  surface  estimates  meander  wavelength  of  approximately  600  km  and 
from  drifters).  The  largest  surface  EKE  values  are  in  maximum  amplitude  of  about  150  km.  Halliwell  and 
the  recirculation  zone  where  the  mean  flow  is  south-  Mooers  (  1983)  found  a  standing  meander  ofcompa- 
westward.  Deep  EKE  is  also  much  weaker  than  ob-  rable  wavelength  downstream  from  Cape  Hatteras  us- 
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ing  Naval  Oceanographic  Office  analyses  of  satellite 
infrared  imagery  from  1975  to  1978.  As  the  Stream 
approaches  the  outflow  boundary  the  meander  ampli¬ 
tude  must  decrease  to  almost  zero.  The  deep  mean 
flow  is  depicted  via  the  mean  pressure  (Fig.  5c),  with 
stationary  closed  circulations  in  phase  with  the  standing 
w  aves  in  the  mean  Gulf  Stream.  There  is  a  strong  baro- 
tropic  component  to  the  circulation  in  these  regions. 
Recirculation  zones  to  the  north  and  south  of  the 
Stream  are  clearly  depicted. 

Experiment  3  is  identical  to  Experiment  2  except 
that  the  entire  eastern  boundary  is  opened.  Both  the 
instantaneous  (Fig.  6a)  and  th**  mean  (Fig.  6b)  path 
of  the  model  Gulf  Stream  are  well  north  of  those  from 
the  narrow  outflow  experiment.  The  amplitude  of  the 
(standing)  meanders  of  the  mean  Stream  in  Experi¬ 
ment  3  are  smaller  in  the  eastern  half  of  the  domain 
than  foi  Experiment  2.  Likewise,  below  the  thermocline 
the  standing  eddies  are  largely  confined  to  the  western 
part  of  the  domain  (Fig.  6c). 

The  tendency  for  the  model  Gulf  Stream  to  extend 
too  far  north  is  very  clear.  One  explanation  for  this 
behavior  is  that  the  Stream  is  not  sufficiently  inertial 
to  separate  at  Cape  Hatteras.  We  have  conducted  sev¬ 
eral  experiments  in  which  we  have  increased  the  max¬ 
imum  current  speed  in  the  Stream  by  up  to  50%,  both 
by  increasing  upper  layer  transport  and  by  using  a 
thinner  upper  layer.  The  Stream  did  not  separate  sig¬ 
nificantly  farther  south  in  these  experiments.  A  large 
barotropic  anticyclonic  recirculation  south  of  the  Gulf 
Stream  forms  just  downstream  from  Cape  Hatteras  in 
the  mean  flow,  tending  to  help  deflect  the  current  north 
along  the  coast.  This  recirculation  intensifies  with  in¬ 
creased  Gulf  Stream  currents  at  inflow.  Surprisingly, 
moving  the  southern  boundary  five  degrees  farther 
south  did  not  significantly  alter  the  tendency  for  over¬ 
shoot  to  occur.  Even  when  the  lateral  boundary  was 
chosen  be  the  actual  coastline  rather  than  the  200  m 
isobath,  overshoot  still  occurred.  While  higher  vertical 
resolution  may  resolve  this  problem,  to  our  knowledge 
all  “realistic"  Gulf  Stream  models  with  higher  vertical 
resolution  share  this  tendency  for  overshoot.  For  ex¬ 
ample  it  is  quite  clear  in  the  results  of  Cox  (  1985.  Fig. 
3b)  in  an  18-level,  one-third  degree  integration  of  the 
GFDl.  model.  As  an  aside.  W.  Holland  (personal 
commur'-  ation  )  finds  that  overshoot  can  be  less  severe 
with  lower  values  of  lateral  eddy  viscosity. 

The  deep  eddy  kinetic  energy  fields  for  experiments 
2  and  3  are  shown  in  Fig.  7a  and  7b.  The  eddy  kinetic 
energy  levels  are  consistently  lower  in  the  narrow  out¬ 
flow  experiment.  Values  near  60  cm2  s  :  occur  in  the 
wide  outflow  case  but  values  less  than  half  that  are 
found  in  the  narrow  outflow  case.  These  two  experi¬ 
ments  present  us  with  a  dilemma.  On  one  hand,  the 
mean  path  is  overall  better  represented  when  we  fix 
the  outflow  boundary  as  a  narrow  port,  but  the  mean 
path  is  too  far  north  when  we  open  the  entire  domain. 


However,  the  deep  eddy  kinetic  energy  is  more  realistic 
in  the  wide  outflow  experiment.  How  might  we  resolve 
these  dichotomous  results? 

Before  answering  this  question  we  should  inquire  as 
to  the  influence  of  topography  on  the  circulation.  Ex¬ 
periment  4  is  defined  by  the  addition  of  the  bottom 
topography  to  Experiment  3.  Following  spinup  the 
model  yields  the  mean  SSH  show  n  in  Fig.  8a.  The  GS 
mean  path  is  quite  unrealistic,  hugging  the  coast  much 
too  far  north  before  turning  eastward.  The  deep  mean 
pressure  in  Fig.  8b  exhibits  length  scales  comparable 
to  those  of  the  topography.  The  mean  anticy  clonic  cir¬ 
culation  around  Bermuda  and  southwest  of  the  south¬ 
ern  portion  of  the  NESC  are  particularly  striking,  as  is 
the  mean  anticyclonic  and  two  cyclonic  gyres  east  of 
the  NESC  underneath  the  Gulf  Stream.  The  deep  EKE 
is  also  very  different  from  that  from  the  flat  bottom 
experiments  (see  Fig.  12a).  With  the  addition  of  bot¬ 
tom  topography  the  maximum  deep  EKE  increases  to 
over  80  cm2  s  1  and  is  largely  confined  west  of  the 
NESC. 

One  approach  to  the  problem  of  separation  of  the 
Gulf  Stream  might  be  to  include  an  additional  vorticity 
source,  such  as  a  deep  southward-flowing  current.  Such 
a  current  was  postulated  by  Stommel  ( 1 958 )  as  a  means 
of  supplying  cold  water  from  the  north  to  feed  the  slow 
upwelling  in  the  ocean  interior  and  was  first  observed 
by  Swallow  and  Worthington  ( 1961 ).  It  has  since  been 
termed  the  Deep  Western  Boundary  Current  (Hogg 
1983).  Estimates  of  its  strength  have  ranged  from  a 
few  Sverdrups  up  to  40  Sv  (Riser  et  al.  1978).  Lai 
(  1 984 )  found  that  an  L-shaped  array  of  current  meters 
along  and  across  the  Blake  Escarpment  yielded  a 
DWBC  with  a  mean  transport  of  24  Sv  over  a  hundred 
day  record.  The  core  of  the  current  was  situated  at 
about  2500  m,  10  km  east  of  the  break  in  the  Escarp¬ 
ment.  Hogg  et  al.  (  1986)  report  a  transport  of  20  Sv 
carried  by  the  DWBC  and  the  slope  current  at  70°W. 
Additionally.  Mellor  et  al.  (  1982 )  found  in  a  diagnostic 
calculation  of  the  general  circulation  of  the  Atlantic  a 
strong  (25  Sx  ).  southwestward  flowing,  nearly  baro¬ 
tropic  current  along  the  continental  slope  in  the  North 
and  Middle  Atlantic  Bight.  Others  et  al.  (  1 985  ).  using 
a  beta-spiral  technique,  have  inferred  a  comparable 
southward  flow  at  2000  m  along  the  continental  slope 
of  the  northwestern  Atlantic.  Mellor  ( personal  com¬ 
munication  )  has  found  this  slope  current  to  be  critical 
for  separation  of  the  Gulf  Stream  oil' Cape  Hatteras  in 
a  sigma-coordinate  primitive  equation  model. 

In  the  model  we  have  imposed  a  DWBC'  as  a  deep 
inflow  along  the  eastern  side  of  the  northern  boundary 
x  ta  a  100  km  wide  stream  centered  near  the  3000  m 
isobath  with  a  parabolic  transport  distribution.  Figure 
9  shows  the  spinup  process  via  the  lime  series  of  po¬ 
tent  al  and  kinetic  energies  averaged  oxer  the  domain. 
The  scale  for  potential  energy  is  fivefold  that  for  kinetic 
energy.  Note  the  time  scale  required  to  reach  statistical 
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Fig.  7.  (a)  EKE2  for  Experiment  2.  (b)  EK.E2  for  Experiment  3.  Contour  interval  is  5  cm2  s'2. 


equilibrium  is  8-10  years  and  that  significant  energetic 
fluctuations  have  time  scales  of  several  months  to  sev¬ 
eral  years  after  equilibrium  is  attained. 

Figure  10a  shows  the  mean  sea  surface  for  Expcri- 
ment  5  (DWBC  =  20  Sv)  after  statistical  equilibrium 
is  reached.  Notice  that  the  mean  Stream  has  separated 
from  the  coast  at  a  lower  latitude  and  that  the  return 
circulation  fills  the  region  to  the  south  of  the  Stream 
in  the  model  interior. 

The  deep  mean  pressure  and  the  deep  mean  velocity 
vectors  are  shown  in  Fig.  10b.  The  mean  velocities 
readily  indicate  the  location  of  the  DWBC  along  the 
continental  slope  and  tend  to  follow  the  contours  of 
f/h.  The  deep  pressure  pattern  is  somewhat  similar  to 
that  without  the  DWBC  but  is  considerably  strength¬ 
ened.  With  the  exception  of  eddy  energy  estimates,  the 
most  extensive  data  set  which  can  be  compared  in  the 


western  North  Atlantic  over  long  time  periods  is  the 
deep  mean  flow.  Hogg  (1983)  compiled  all  of  the  long¬ 
term  deep  current  meter  measurements  in  this  region. 
He  has  recently  updated  that  map  (Hogg  et  al.  1986) 
as  shown  in  Fig.  lc.  While  there  are  numerous  detailed 
comparisons  to  be  made  we  w  ill  focus  on  four  specific 
areas:  1 )  the  recirculation  /one  near  the  region  of  the 
HEBBLE  experiment  at  41°N,  63° W.  2)  the  foot  of 
the  continental  rise  near  68°W,  37°N.  3)  the  region 
of  the  “northern  recirculating  gyre"  north  of  40°N 
along  55°W.  and  4 )  the  flow  around  isolated  seamounts 
such  as  Corner  Rise. 

Hogg  ( 1983)  identified  two  separate  deep  flow  re¬ 
gimes,  the  classical  DWBC  flowing  along  the  conti¬ 
nental  rise  and  a  second  system  of  recirculating  classical 
slope  water  gyres  farther  offshore.  However  these  re¬ 
gimes  are  not  always  distinct,  for  example,  Hogg 
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4800  m  isobath  near  37°N,  68°W  is  found  in  the  ob¬ 
servations.  The  model  results  clearly  point  toward  the 
need  for  additional  measurements  to  determine  if  a 
closed  gyre  exists. 

In  Hogg’s  1983  schematic  the  current  meter  array 
near  40°N,  59°  W  was  not  available  and  Hogg  suggested 
the  deep  flow  in  this  region  was  westward.  However. 
Fig.  1c  (from  1986)  clearly  shows  that  the  currents  in 
this  area  are  flowing  toward  the  northeast.  The  model 
results  show  that  the  eastern  side  of  the  “Hebble  recir¬ 
culation"  is  also  the  western  side  of  a  small  antieydonic 
gyre  which  splits  what  Hogg  termed  the  “northern  re¬ 
circulating  gyre"  into  two  separate  closed  cyclonic 
gyres.  The  anticyclonic  circulation  is  reminiscent  of 
that  suggested  by  Schmitz  ( 1980)  [and  by  Richardson 
(1981)  from  one  degree  box -averaged  temperature  data 


( 1983)  and  Hogg  et  al.  ( 1986)  found  that  the  Wunsch 
and  Grant  ( 1982)  circulation  scheme,  deduced  from 
an  inversion  of  hydrographic  data,  could  not  easily  ac¬ 
count  for  the  HEBBLE  data.  It  suggested  as  much  as 
40  Sv  flowing  along  the  4000  m  isobath  near  63°W 
was  not  evident  at  70°W.  Hogg  surmised  that  there 
was  a  decoupling  between  the  two  currents  such  that 
a  recirculation  existed  in  the  deep  flow  in  the  HEBBLE 
area.  Such  a  recirculation  is  very  much  evident  in  the 
deep  mean  flow  from  the  model  with  the  magnitude 
of  the  recirculation  being  about  35  Sv.  This  circulation 
is  indeed  decoupled  from  that  near  70°W.  In  fact,  a 
separate  “near-mesoscale"  recirculation  zone  exists  to 
the  southwest  along  the  foot  of  the  continental  rise 
centered  near  37°N.  68°W  and  may  explain  the  Rise 
Array  data.  In  particular,  the  shoreward  flow  at  the 
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Fig  .  9.  Time  series  of  basin-averaged  potential  energy  per  unit 
area  due  to  thermocline  (model  interface)  fluctuations  and  upper 
and  lower  layer  kinetic  energy  per  unit  area  in  units  of  kg  s~J.  Scale 
on  the  left  refers  to  kinetic  energy  while  scale  on  the  right  refers  to 
potential  energy. 


at  450  m],  although  the  model  gyre  is  farther  north 
and  west  of  that  indicated  in  the  observations.  The 
eastern  cyclonic  gyre  centered  near  41  °N,  55°  W  is  dy¬ 
namically  similar  to  that  outlined  by  Hogg  and  Stom- 
mel  ( 1985)  in  that  this  deep  gyre  is  nearly  potential 
vorticity  conserving. 

To  illustrate  the  influence  of  the  DWBC  on  the  Gulf 
Stream  we  have  plotted  in  Fig.  1 1  the  observed  mean 
path  of  the  15°C  isotherm  at  200  m  (adapted  from 
Fisher  1977 )  versus  the  axis  of  the  mean  current  above 
the  thermocline  for  three  experiments  differing  only  in 
the  intensity  of  the  lower  layer  flow:  No  DWBC,  a  20 
Sv  DWBC,  and  a  40  Sv  DWBC.  As  we  increase  the 
DWBC  magnitude  the  mean  path  of  the  model  jet  is 
shifted  further  south.  In  fact  the  separation  latitude  is 
shifted  by  1 50  km  over  the  range  of  transports.  As  the 
transport  of  the  DWBC  is  increased,  the  model  Gulf 
Stream  path  is  shifted  southward,  agreeing  best  with 
observations  at  a  transport  of  20  Sv.  West  of  the  sea¬ 
mounts,  as  the  DWBC  intensifies,  the  standing  mean¬ 
der  amplitude  and  the  recirculation  transport  tend  to 
decrease.  Near  Cape  Hatteras,  there  is  an  overshoot  of 
all  modeled  paths  relative  to  observation,  decreasing 
in  amplitude  and  zonal  scale  with  increasing  transport. 
Note  that  the  mean  current  is  mainly  north  of  the 
NESC  in  Expts.  4  and  5  but  through  them  in  Experi¬ 
ment  6.  Interestingly,  Auer  ( 1987)  has  found  a  statis¬ 
tically  significant  northward  shift  in  the  Gulf  Stream 
over  five  years  from  1980  to  1985  of  more  than  half  a 
degree  latitude.  He  has  no  explanation  for  the  shift. 

What  are  the  dynamics  which  allow  the  DWBC  to 
influence  the  mean  path  of  the  Gulf  Stream?  For  the 
two-layer  system  they  are  very  similar  to  those  discussed 
in  HT80  to  explain  how  deep  southward  flow  along 


the  west  Florida  shelf  could  prevent  northward  pene¬ 
tration  of  the  Loop  current.  HT84  also  used  these  dy¬ 
namics  to  show  how  fluctuations  in  the  deep  flow,  rec¬ 
tified  by  the  NESC,  could  alter  the  mean  path  of  the 
Gulf  Stream.  From  simple  kinematics,  the  continuity 
equation  for  the  upper  layer  can  be  written  as 

Ai,  + /tiV  •  v,  +  v,  •  VA|  =  0.  (3.1) 

The  third  term  on  the  left  hand  side  of  ( 3. 1 )  represents 
nonlinear  advection  of  upper  layer  thickness  by  a  com¬ 
ponent  of  the  upper  layer  flow.  For  nearly  geostrophic 
flow, 

Vi-VA,  «=  v,g*  VA,  =  v2g-Vh,  =  J(th,  r?2),  (3.2) 

where  vIJ?  is  the  geostrophic  velocity  component  in  layer 
/,  J  is  the  Jacobian  operator,  r/,  is  the  deviation  of  the 
free  surface  from  its  initial  value,  and  ?/2  is  the  deviation 
of  the  interface  between  layers  from  its  initial  value. 
The  term  can  be  large  when  a  strongly  tilted  thermo¬ 
cline  (the  interface  in  the  Gulf  Stream)  intersects  an 
intense  deep  flow  (the  DWBC)  at  a  large  angle.  Phys¬ 
ically  it  represents  the  interaction  of  a  barotropic  cur¬ 
rent  and  baroclinic  current.  Thus  the  barotropic  com¬ 
ponent  of  the  DWBC  can  advect  the  Gulf  Stream  at 
the  point  of  intersection  near  the  continental  slope. 
For  example,  in  Fig.  1 1  let  the  0  DWBC  curve  represent 
the  steady  state  position  of  the  Gulf  Stream.  If  we  in¬ 
crease  the  DWBC  transport  to  20  Sv.  say,  nearly  geo¬ 
strophic  deep  flow  moves  southward  closely  following 
//A  contours  (dominated  by  bottom  topography  nearly 
everywhere).  Eventually  the  DWBC  flows  underneath 
the  Gulf  Stream.  In  Eq.  (3.1 )  the  initial  steady  state 
balance  (with  A„  =  0)  is  altered  such  that,  in  the  region 
of  intersection  of  the  Gulf  Stream  and  the  DWBC. 
v,  •  Vfi,  =  v2lth\y>  0,  since  v2g  <  0  and  hiy  <  0.  Thus 
A„  <  0  and  the  upper  layer  thickness  decreases  locally. 
As  the  layer  thins,  potential  vorticity  is  conserved  via 
the  flow  gaining  anticyclonic  relative  vorticity  and/or 
leaving  the  coast  at  a  more  southerly  latitude,  in  effect 
shifting  the  Gulf  Stream  southward  along  the  conti¬ 
nental  slope,  as  indicated  in  the  20  Sv  (dashed-dot) 
curve  of  Fig.  1 1 . 

The  significant  influence  of  the  DWBC  on  the  Gulf 
Stream  path  clearly  extends  into  the  model  interior. 
This  change  in  the  interior  has  radical  consequences 
for  the  deep  eddy  field  as  well.  In  Fig.  1 2  we  have  plotted 
the  deep  eddy  kinetic  energy  for  the  0-,  20-.  40-Sv 
DWBC  experiments  (4.  5a,  6  respectively).  In  the  case 
with  no  DWBC  (Fig.  1 2a)  the  deep  EKE  is  largely 
confined  to  a  region  west  of  the  NESC  with  maximum 
values  less  than  80  cm2  s  :.  The  NESC  is  also  evident 
in  the  EKE  contour  pattern.  In  Fig.  1 2b  we  have  plotted 
the  EKE  for  the  20-Sv  experiment.  The  distribution  of 
EKE  is  more  equally  distributed  east  and  west  of  the 
NESC  with  maximum  values  near  50  enr  s  '  to  the 
east  and  90  cm 2  s " 2  to  the  west.  For  the  40-Sv  DWBC 
case  the  EKE  field  is  almost  entirely  shifted  to  east  of 
the  NESC  with  maximum  values  near  100  cm2  s''2. 
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Fig.  10.  As  in  Fig.  8  but  for  Experiment  5a. 


These  three  experiments  demonstrate  that  increasing 
the  DWBC  shifts  the  EKE  distribution  further  eastward 
and  toward  higher  amplitude  as  the  mean  path  of  the 
Gulf  Stream  shifts  southward. 

As  we  have  noted  earlier,  an  extensive  comparison 
of  long-term  statistics  for  models  and  observations  was 
conducted  by  HS86,  primarily  for  deep  EKE  in  the 
western  North  Atlantic  and  Pacific.  This  type  of  model/ 
data  intercomparison  will  continue  to  be  a  key  test  for 
any  ocean  model  designed  to  simulate  and  explain  the 
dynamics  of  midlatitude  western  boundary  currents 
and  eddies.  Figure  lb  shows  a  new  updated  map  of 
deep  eddy  kinetic  energy  in  the  western  North  Atlantic. 
This  is  to  be  compared  with  the  contour  maps  of  deep 
EKE  provided  for  Experiments  4,  5  and  6  in  Fig.  12a- 
c,  respectively.  As  mentioned  earlier,  as  the  DWBC 


transport  increases  the  area  of  high  deep  EKE  shifts  to 
the  east.  For  no  DWBC  the  maximum  value  of  about 
85  cm2  s  2  is  located  near  37°N.  67°W,  west  of  the 
seamounts,  and  is  comparable  to  the  observed  values. 
Note  also  that  large  values  also  occur  near  the  northern 
end  of  the  NESC  but  confined  to  a  very  small  area. 
Some  indication  of  this  same  small  scale  distribution 
is  shown  in  Fig.  lb  as  well,  with  values  of  1 19  cm2  s  2 
observed  in  a  comparable  location  to  that  in  the  model. 
Near  the  model  boundary  there  are  also  narrow  max¬ 
ima  in  EKE  due  to  topographic  waxes.  East  of  the 
NESC  there  are  no  values  in  excess  of  40  cm  s  2  for 
Experiment  4.  When  the  DWBC  is  increased  to  20  Sv 
in  Experiment  6  the  deep  EKE  field  exhibits  maxima 
to  the  east  and  west  of  the  seamounts,  with  those  to 
the  west  near  68.5°W  being  as  large  as  90  cm2  s'2  and 
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Fig.  11.  Paths  of  the  mean  Gulf  Stream.  Axes  of  the  upper  layer  current  maximum  for  the 
long  term  mean  for  Experiments  4,  Sa,  and  6  and  an  estimate  of  the  mean  path  of  the  Gulf 
Stream  as  indicated  by  the  mean  location  of  the  I5°C.  Isotherm  at  200  m  from  Fisher  (  1977). 


those  to  the  east  near  55°W  not  exceeding  50  cm2  s'2. 
Finally,  for  Experiment  6  with  40  Sv  for  the  DWBC 
the  eastward  shift  of  the  EKE  field  is  nearly  complete, 
having  maximum  values  greater  than  100  cm2  s'2  near 
39°N,  53°W.  This  intercomparison  of  deep  EKE  sug¬ 
gests  that  the  model  amplitudes  are  too  small  by 
roughly  a  factor  of  two  in  comparison  with  the  obser¬ 
vations.  All  else  being  the  same,  time-dependent  fluc¬ 
tuations  in  the  inflow  transport  (or  angle)  and  the  wind 
driving  should  both  act  to  increase  deep  EKE.  While 
fluctuations  in  direct  wind  driving  are  unlikely  to  con¬ 
tribute  significantly  to  abyssal  EKE  (Niiler  and  Kob- 
linsky  1985),  fluctuations  in  volume  transport  on  sea¬ 
sonal  to  interannual  time  scales  near  Cape  Hatteras 
may  be  much  more  critical.  We  do  not  suggest  that  the 
vertical  and  horizontal  resolution,  friction  parameters, 
and  stratification  are  less  important  in  altering  the  am¬ 
plitude  and  distribution  of  EKE  in  the  model.  In  fact 
model  sensitivity  studies  suggest  that  they  may  be  more 
important  than  time-dependent  forcing  in  this  regard. 

In  Fig.  13a  we  show  the  rms  SSH  variability  of  Ex¬ 
periment  5a  (20-Sv  DWBC).  The  maximum  value  is 
about  40  cm  to  the  west  of  the  NESC  chain.  This  cor¬ 
responds  well  in  amplitude  with  values  determined  by 
Marsh  et  al.  ( 1984)  based  on  SEASAT  and  GEOS-3 
crossover  data  although  location  of  the  maximum  is 
several  degrees  too  far  west  in  the  model  results.  With 
the  advent  of  GEOSAT  we  will  have  additional  data 
on  variability  of  the  SSH  in  the  Gulf  Stream  region  to 
compare  with  these  model  results. 

Figure  1 3b  shows  the  associated  field  of  thermocline 
variability  (rms)  for  Experiment  5a.  The  magnitudes 
are  comparable  to  those  observed.  For  example.  Fig. 
13c  compares  the  maximum  rms  variability  from  Ex¬ 
periment  5b  at  each  longitude  with  Dantzler’s  ( 1977) 
calculation  of  the  rms  variability  of  the  1 5°C  isotherm 
at  37°  and  39°N.  The  magnitude  of  the  rms  thermo¬ 


cline  variability  is  10-20%  higher  in  the  experiment 
with  bottom  friction  (5b)  than  without  (5a).  We  are 
encouraged  by  two  aspects  of  this  comparison.  First, 
the  amplitudes  of  rms  thermocline  variability  of  the 
model  are  only  slightly  lower  over  the  range  of  longi¬ 
tudes  than  observed.  Since  our  model  thermocline  has 
a  mean  depth  of  1000  m  the  variability  of  the  15°C 
isotherm  is  probably  underestimated.  This  is  offset  by 
the  fact  that  model  variability  was  chosen  not  along 
constant  latitudes  lines  but  the  latitude  of  maximum 
variability.  Second,  the  model  rms  variability  did  not 
decrease  significantly  downstream  to  45°N.  at  the  out¬ 
flow  boundary.  This  suggests  the  boundary  condition 
did  not  artificially  constrain  stream  variability  near  the 
outflow. 

The  shallow  EKE  estimate  from  drifter  data  is  shown 
in  contour  form  in  Fig.  la  (Schmitz  et  al.  1983:  Rich¬ 
ardson  1981).  Model  values  are  provided  in  Fig.  14 
for  the  upper-layer  EKE  for  Experiments  5a  (a)  and 
5b  (b).  Note  the  upper-layer  EKE  levels  to  the  east  of 
the  NESC  are  higher  for  the  experiment  with  bottom 
friction  although  not  very  different  to  the  west  of  the 
seamounts.  A  serious  consideration  with  this  compar¬ 
ison  is  that  the  mean  thickness  of  the  model  upper 
layer  is  1000  m  while  the  surface  drifter  measurements 
are  representative  of  the  upper  100  meters  or  so.  Based 
on  estimates  of  vertical  mean  shear  in  the  Gulf  Stream 
region  the  model  values  are  reasonable  for  the  vertical 
average  over  1000  m  and  are  roughly  half  the  surface 
EKE  values  (see  SH86).  The  maximum  value  from 
observations  is  about  3000  cm2  s'2  located  near  67° W 
while  in  the  model  the  maximum  value  is  about  2400 
centered  near  the  same  longitude.  In  Fig.  14b  the  area 
covered  by  variability  greater  than  1000  cm2  s  2  is 
considerably  greater  than  that  from  the  model  results. 
The  downstream  extent  of  the  high  variability  in  the 
model  is  comparable  to  that  from  the  observations. 
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although  l Ik*  downstream  decay  appears  more  rapid  in 
the  drifter  data  than  in  the  model  results.  As  discussed 
above,  the  deep  layer  thickness,  the  lar  k  ol  time  de¬ 
pendence  in  the  wind  field  and  the  constant  inflow 
transport  all  contribute  to  the  lower  values  in  I  KE 
from  the  model  results.  It  is  very  likely  that  highei  \cr 
lical  and  horizontal  resolution  in  the  model  would  im¬ 
prove  the  comparison  with  the  data. 

f  igure  15  shows  a  histogram  for  the  mass  transport 
at  55" W  for  the  (a)  upper  layer.  ( b )  lower  layer,  and 
( c )  total,  averaged  over  one  degree  latitude  bands.  I  he 
recirculation  zones  are  dramatically  illustrated  in  these 
figures.  Note  that  the  total  eastward  transport,  which 
includes  a  significant  eddy-flux  induced  circulation,  is 
194  Sv  between  59°  and  41  ”N  ( recall  the  specified  in¬ 
flow  at  Cape  l  ear  is  50  Sv).  This  is  much  larger  than 
the  total  mean  transport  of  95  Sv  at  55°W  by  Rich¬ 


ardson  (  1985).  However,  the  appropriately  weighted 
difference  between  model  upper  and  lower  layer  trans¬ 
ports  associated  with  the  model  (  mil  Stream  analogue 
v  ielils  a  value  ( relative  to  the  bottom  )  of  95  Sv .  close 
to  the  90  Sv  estimate  In  McCartney  et  al  (  19X0)  for 
the  same  calculation  from  observations  at  55''W.  (See 
also  Schmitz  1980). 

4.  Summary  and  conclusions 

We  have  applied  a  prim  it  iv  c-equation  numerical 
model  as  a  two-layer  analogue  of  the  ( iulf  Stieam  Sys¬ 
tem  to  a  limited  area  from  (  ape  I  latteras  to  east  of  the 
(iraiul  Banks  (  78"-45"W,  50"  4X"N).  In  this  initial 
paper  six  numerical  experiments  from  a  series  of  over 
50  integrated  to  statistical  ctpiihhrium  were  selected 
for  detailed  description  and  intercomparison  with  ob- 
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Fro.  15.  Mean  zonal  transport  histogram  averaged  over  2  degrees  latitude  hands  at  55°W.  from 
Experiment  5a  for  the  (a)  upper  layer,  (b)  lower  layer,  and  (c)  total.  Each  degree  longitude  represents 
10  Sv  of  transport. 
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servntion.  Experiment  1  consisted  of  a  flat  bottom  re-  scribed  for  two  different  outflow  specifications  at  the 
gime  driven  by  wind  forcing  only.  Realistic  inflow  eastern  boundary  in  Experiments  2  and  3.  Three  other 
transport  above  the  main  thermoclinc  was  then  pre-  model  runs  included  (4)  bottom  topography.  (5)  a  deep 
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western  boundary  current  with  20  Sv  total  transport 
added  to  Experiment  4,  and  (6)  increasing  the  Deep 
Western  Boundary  Current  ( DWBC )  in  Experiment  5 
to  40  Sv. 

Without  a  prescribed  inflow  (Florida  Current 
Transport)  the  resulting  velocity  fields  are  much  too 
weak,  as  expected.  For  a  realistic  inflow  transport  ex¬ 
iting  at  a  fixed  narrow  outflow  location  both  the  mean 
and  eddy  fields  are  quite  unrealistic.  For  an  entirely 
open  eastern  boundary  the  modelled  Gulf  Stream  Sys¬ 
tem  with  realistic  topography  is  more  energetic  at 
abyssal  depths  but  is  located  too  far  north  relative  to 
observations.  The  addition  of  the  DWBC  moves  the 
mean  path  of  the  midlatitude  jet  southward  toward  a 
more  realistic  configuration,  shifting  the  region  of  in¬ 
creased  abyssal  EKE  eastward  thus  improving  the 
comparison  with  deep  observations  at  55°W.  The  20 
Sv  DWBC  experiment  (Expt.  5)  also  yields  a  deep 
mean  flow  which  corresponds  to  the  observational  data 
in  several  areas,  including  recirculations  north  and 
south  of  the  Gulf  Stream  and  a  deep  cyclonic  gyre  just 
to  the  east  of  the  northern  part  of  the  NESC.  Another 
deep  cyclonic  circulation  at  the  foot  of  the  continental 
rise  centered  near  37°N,  68°W  and  anticyclonic  deep 
flows  around  several  isolated  topographic  features  such 
as  the  Corner  Rise  Seamounts  and  Bermuda  are  also 
indicated  in  the  model  results  but  as  yet  have  little  or 
no  observational  data  to  confirm  them.  The  amplitudes 
of  upper  level  eddy  kinetic  energies  and  the  mean  flow 
are  within  at  least  a  factor  of  two  of  the  data  base  and 
the  model  equivalent  Gulf  Stream  transport  at  55°W 
is  potentially  relevant.  Additional  model/data  inter¬ 
comparisons  utilizing  GEOSAT  altimeter  data,  in¬ 
verted  echo  sounders,  and  bottom  pressure  gauges  are 
reported  by  Hallock  et  al.  ( 1988). 

These  initial  results  are  encouraging,  but  the  need 
for  future  work  is  clear.  To  be  examined  further  are 
problems  associated  with  changing  vertical  and  hori¬ 
zontal  resolution  (Are  our  two-layer  results  valid  in 
the  limit  of  continuous  stratification?),  model  dynam¬ 
ics  ( How  is  the  Northern  Recirculation  Gyre  driven 
by  eddy  fluxes?),  and  model/data  intercomparisons 
[Are  model  simulations  of  warm  and  cold  ring  gen¬ 
eration.  translation,  and  decay  compatible  with  the 
observations  of  Auer  (  1987 ).  Brown  el  al.  (  1986).  and 
others?].  We  do  not  discount  the  possibility  that  sep¬ 
aration  of  the  Gulf  Stream  could  be  better  achieved  by 
allowing  surfacing  of  the  fluid  interface.  A  surfaced  or 
ventilated  thermoclinc  appears  to  be  a  crucial  com¬ 
ponent  in  a  successful  thermoclinc  theory  ( Luyten  et 
al.  1983),  and  Veronis  (1981)  have  argued  that  realistic 
separation  of  midlatitude  western  boundary  currents 
from  the  coast  must  necessarily  include  surfacing  of 
the  thermoclinc  as  well.  Western  boundary  currents 
forced  out  to  sea  between  wind  gyres  of  opposite  sign 
are,  in  this  view,  qualitatively  different  because  the 
thermocline  does  not  surface.  Nevertheless,  we  believe 
the  key  issue  discussed  in  this  paper,  the  influence  of 


the  DWBC  on  the  Gulf  Stream,  is  relevant  to  either 
regime.  Because  of  the  efficiency  and  relative  simplicity 
of  this  model  we  have  begun  to  examine  these  and 
related  problems  using  computational  resources  widely 
available  today. 
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APPENDIX  A 

List  of  Symbols 

A  horizontal  eddy  viscosity 

C),  coefficient  of  bottom  friction 

/  Coriolis  parameter 

£  acceleration  due  to  gravity 

g'  reduced  gravity,  equal  to  g(  i>i  ~  t>\)/ 1> 

II,  thickness  of  layer  i  at  rest 

h,  instantaneous  layer  thickness  of  the  /  th  layer 

p,  density  normalized  pressure  for  layer  / 

t  time 

v,  current  velocity  in  layer  /  with  v-directed 

component  it,  and  r-directed  component 

v, 

V,  volume  transport.  h,v, 

\,K  geostrophic  current  velocity 

v,  y,  r  tangent  plane  Cartesian  coordinates:  v  pos¬ 
itive  eastward,  y  positive  northward,  and 
r  positive  upward 

A/  time  step  in  the  numerical  integration 
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Ax ,  Ay  horizontal  grid  increments  for  each  inde¬ 
pendent  variable 

ih  free  surface  height  anomaly  (FSA),  height 

of  the  free  surface  above  its  initial  eleva¬ 
tion 

p,  pi  densities  of  seawater 
t,x,  Tjy  x-  and  >-directed  tangential  stresses  at  the 
top  (/)  and  bottom  (/'  4  1 )  of  layer  i  re¬ 
spectively 


APPENDIX  B 

The  Gulf  Stream  Model 


Consider  a  stably  stratified  Boussinesq  fluid  with  a 
free  surface  and  a  fixed  density  contrast  specified  be¬ 
tween  immiscible  layers.  All  thermodynamic  effects  are 
neglected.  Each  interface  between  fluid  layers  is  thus 
an  isopycnal  surface.  While  the  actual  system  of  equa¬ 
tions  is  solved  on  an  earth-sized  spherical  geometry, 
for  simplicity  we  will  use  Cartesian  geometry  in  the 
discussion  below.  In  transport  form  the  vertically  in¬ 
tegrated  equations  of  motion  used  in  the  models  are 

dV 

~  4  (V-V,  4  V, •  V ) v,  4  k  X  /  V, 
at’ 


=  -h,VP,  +  p  '(r,  r,.,)  4  .-tC-’V,.  (B.l) 


where 


at 


4 


V 


V,  =  0. 


V 


a_ 

dv 


j. 


/>, 


l  I 

=  ft?  Vi  H  Z 

i  l 


4/  -\P,~  p,)/p,. 

v,  =  //,v,  =  h a  It, i  4  V,i), 

r,  =  r,'i  4  t,'  j. 


( B.2 ) 


( B.3 ) 


See  appendix  A  for  symbol  definitions. 

In  the  two-layer  experiment  discussed  here,  i  =  1.2 
with  t2  (interfacial  stress)  set  to  zero  and  t,  (bottom 
stress)  quadratic  in  velocity,  i.e. 


rj  C,  !  v-.  I  v:. 

This  layered  model,  with  transports  L\  I  ' as  depen¬ 
dent  variables,  handles  strongly  sloping  topography  es¬ 
pecially  well,  so  long  as  it  is  confined  to  the  lowest 
layer.  The  topography  appears  multiplicatively  in  the 
pressure  gradient  term  and  is  differentiated  only  to  the 
extent  that  it  affects  the  velocity  field  in  the  advective 
terms.  When  large-amplitude  topography  is  introduced, 
restrictions  on  the  time  step  and  the  eddy  coefficient 
are  affected  only  to  the  extent  that  the  topography  de¬ 


termines  the  amplitude  of  the  velocity  field  (advective 
CFL  condition)  and  the  grid  interval  Reynolds  number. 
The  continuity  equation  is  linear  when  written  in 
transport  form,  thus  avoiding  complications  from  the 
nonlinear  advective  term  when  layer  thickness  and  ve¬ 
locity  are  used  as  variables. 

Except  at  the  inflow  and  outflow  ports,  the  walls  are 
rigid  and  the  no-slip  condition  is  prescribed  on  the 
tangential  flow.  At  inflow  the  profile  of  transport  V,  is 
prescribed  as  a  parabolic  function.  Because  the  current 
is  nearly  geostrophic,  the  upper  layer  velocity  maxi¬ 
mum  is  to  the  left  of  the  center  of  the  jet  when  facing 
downstream.  The  normal  flow  at  the  outflow  port  is 
self-determined  with  the  integral  constraint  that  the 
total  transport  out  from  each  layer  match  the  inflow 
at  all  times. 

The  numerical  design  of  the  model  is  similar  to  that 
described  by  Hurlburt  and  Thompson  ( 1980)  and  will 
not  be  rewritten  here.  However,  differences  in  meth¬ 
odology  for  handling  the  open  boundary  conditions 
are  described  below.  The  flow  at  the  open  boundaries 
are  predicted  in  three  steps: 

First,  a  modified  Orlanski  ( 1976)  radiation  condition 
is  used  such  that  for  each  component  of  the  transport 


( B.4 ) 


4* 


'2  At 

+  I  —  |<  </«  i 
A  A 


Ate 

4  — 

Ax 


where  B  refers  to  the  outflow  boundary  .  B  -  1  is  one 
grid  point  interior  to  the  boundary. 


(  =  min 


Ax 

At 


.  max  (a.  uM ). 


tl,n  is  the  mean  inflow  speed,  and  the  phase  speed  is 
determined  locally  for  each  wave  mode.  There  are  sev¬ 
eral  ways  to  determine  the  local  phase  speed.  <■*.,  but 
we  have  found  that  <<  =  iiM  is  effective  for  outflow  for 
the  Gulf  Stream. 

Second,  if  inflow  is  predicted  at  an  outflow  point,  a 
linear  implicit  drag  is  applied.  This  is  required  because 
in  an  open  boundary  condition  with  self-determined 
flow  the  mean  is  constrained.  In  the  case  of  an  east- 
west  boundary,  deviations  from  the  mean  flow  are 
somewhat  constrained  by  the  tendency  to  form  a  Sv  er¬ 
drup  interior  (Hurlburt  and  Thompson  1973).  How¬ 
ever.  for  a  north-south  oriented  boundary  it  is  possible 
to  develop  a  recirculation  at  the  eastern  open  boundary 
driven  by  the  eddy-fluxes  from  the  Gulf  Stream.  This 
mode  can  grow  unrealistically  large  unless  there  is  a 
“brake"  on  the  inflow  component  of  the  t  circulation 
at  the  outflow  boundary  (Hurlburt  and  Thompson 
1980).  The  damping  time  scale  (one  day)  effectively 
inhibits  an  unrealistic  recirculation  at  45°W. 

Finally,  the  x-component  of  transport  is  uniformly 
accelerated  or  deccclerated  such  that  total  inflow  mass 
transport  is  matched  by  total  outflow  transport  for  that 
time  step.  The  outflow  boundary  condition  docs  not 
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use  any  interior  values  at  time  level  n  +  1,  so  it  can 
provide  the  solver  with  the  boundary  information  on 
transport  required  at  the  ports.  We  can  then  solve  the 
Helmholtz  equations  for  h 1 .  Once  we  have  h") 1 ,  the 
semi-implicit  calculation  of  L'?*\'/2.j  and  I'"]1 ,  /2  in  the 
interior  can  be  completed. 

The  behavior  of  the  open  boundary  condition  has 
been  remarkably  good.  Comparisons  by  Heburn  and 
Hurlburt  (personal  communication )  of  results  from  a 
Caribbean  model  with  the  radiation  condition  used 
here  and  an  identical  experiment  using  the  boundary 
condition  of  Hurlburt  and  Thompson  ( 1980)  without 
the  Coriolis  term  show  good  agreement.  Strong  jets 
pass  through,  as  do  eddies.  Gravity  wave  reflections 
are  not  a  dominant  component  of  the  model  solutions 
nor  is  extra  damping  required  on  the  outflow  boundary. 

The  Helmholtz  equations  are  solved  using  a  vecto¬ 
rized  variant  of  the  capacitance  matrix  technique 
(Busbee  et  al.  1971)  for  homogeneous  Neumann 
boundary  conditions  one-half  grid  distance  outside  the 
h  grid  on  a  general  non-rectangular  geometry  ( Wallcraft 
1980).  Each  solution  requires  two  calls  to  a  variant  of 
Hockney’s  ( 1965 )  method  developed  hv  Daniel  Moore 
(Imperial  College,  London)  and  Alan  Wallcraft  for 
rectangular  regions. 

REFERENCES 

Adamec.  D..  1988:  Numerical  simulations  of  the  effects  of  seamounts 
and  vertical  resolution  on  strong  ocean  flows.  J  Phvs  Oceanogr . 
18,  258-269. 

Auer.  S.  J..  1987:  Five-year  climatological  survey  of  the  Gulf  Stream 
Svstem  and  its  associated  rings  J  (ieophvs  Res  92,  1 1  709- 
11  726. 

Bogue.  N.  M  .  R  X.  Huang  and  K.  Bryan.  1986:  Verification  ex¬ 
periments  with  an  Isopscnal  coordinate  model  J  Pin  s  Rev  . 
16.  985-990. 

Boris.  J.  P..  and  D.  1..  Book.  197V  Flux-corrected  transport  Part  1: 
SHASTA.  A  fluid  transport  algorithm  that  works  J  Com  pul 
Phvs .  II,  38-69 

Brown.  O.  B  .  P  C.  Cornillon.  S.  R  F.mmerson  and  H  M  Carle. 
19X6:  Gulf  Stream  warm  rings:  a  statistical  studs  of  their  be¬ 
havior  Deep-Sea  Res  .  3.3,  1459-1473 
Busbee.  B  I...  F.  W  Doit.  J  A  George  and  G  IF  Golub.  I97|  fhe 
direct  solution  of  the  discrete  Poisson  equation  on  irregular  re¬ 
gions  .57  13/  ./  Sum  Anal.  722-736. 

Cox.  M  D  ,  19X5.  An  eddy  resolving  numerical  model  of  the  ven¬ 
tilated  thermoclinc.  J  Pin  s  Oceanogr .  10,  1312-1324 
Dant/ler.  H.  I  ..  1977.  Potential  energy  maxima  in  the  tropical  and 
subtropical  North  Atlantic.  J  Pins  Oceanogr  7,512-519. 
Douglas.  B  R.  F  Cheney  and  R  W  Agreen.  19X3:  Eddy  energy 
of  the  Northwest  Atlantic  and  Gulf  of  Mexico  determined  from 
GFOS  3  altimetry  J  (ieophvs  Res  .  88.  9595-9604 
Fisher.  A..  1977:  Historical  limits  of  the  northern  edge  of  the  Gull 
Stream  GuUstream.  3.  265-271 

Halliwell.  G  R  .  Jr  .  and  C.  N  k  Moocrs.  19X3  Meanders  of  the 
Gulf  Stream  downstream  from  Cape  I  lattcras  ./  Pins  (Kean 
ngr  13,1275-1292. 

Hallock.  Z  R  .  J  1.  Mitchell  and  J  I)  Thompson.  I9SS  Sea  surlace 
height  vanability  near  the  New  England  Seamounts:  an  inter- 
comparison  among  in  situ  observations.  GFOS  AT  aiunieuy. 
and  numencal  simulation  ./  (ieophvs  Res  in  press. 

Harrison,  D  F...  1982:  On  deep  mean  flow  generation  mechanisms 
and  the  abyssal  circulation  of  numencal  model  gy  res  Dyn  At¬ 
mos  Oceans.  6,  135-152. 


Heburn,  G.  W.,  and  H.  E  Hurlburt.  1988:  A  numencal  study  of 
mesoscale  variability  in  the  eastern  Caribbean  Sea  (submitted 
to  J  Phys.  Oceanogr  ). 

Hellerman.  S..  and  M.  Rosenstein,  1983:  Norma!  monthly  wind  stress 
over  the  world  ocean  with  error  estimates.  J  Phvs  Oceanogr 
13,  1093-1 104. 

Hockney.  R.  W..  1965:  A  fast  direct  solution  of  Poisson's  equation 
using  Fourier  Analysis.  J  Assoc  Comput  Mach  .  12,  95-1 13 
Hogg.  N.  G..  1983:  A  note  on  the  deep  circulation  of  the  western 
North  Atlantic:  Its  nature  and  causes.  Deep-Sea  Res  .  30,  945- 
961. 

- ,  and  H.  Stommel,  1985:  On  the  relation  between  the  deep  cir¬ 
culation  and  the  Gulf  Stream.  Deep-Sea  Res  32,  1181-1193. 

— — ,  R.  S.  Pickart,  R  M  Hendry  and  W.  J.  Smethie.  Jr..  1986.  The 
northern  recirculation  gvre  of  the  Gulf  Stream.  Deep-Sea  Res . 
33,1139-1165. 

Holland.  W.  R.,  and  L.  B.  Fin.  1975:  On  the  generation  of  mesoscale 
eddies  and  their  contribution  to  the  oceanic  general  circulation. 
Part  I:  A  preliminary  numencal  experiment.  J  Phvs  Oceanogr.. 
S,  642-657 

— — ,  and  W.  J.  Schmitz.  1985:  Zonal  penetration  scale  of  midlatitude 
jets.  J  Phys  Oceanogr  15,  1859-1875. 

Hurlburt.  H.  F...  1984:  The  potential  for  ocean  prediction  and  the 
role  of  altimeter  data.  Mar  (Send..  8,  17-66. 

- .  and  J  D.  Thompson.  1980:  A  numerical  study  of  Loop  Current 

intrusions  and  eddy  shedding.  J  Phys  Oceanogr  .  10,  1611- 
1651. 

- ,  and - ,  1984:  Preliminary  results  from  a  numencal  study 

of  the  New  England  Seamount  chain  influence  on  the  Gulf 
Stream.  Pros  of  the  Workshop  on  Predictability  of  Fluid  Motions. 
Amer.  Inst.  Phys.  G.  Hollaway,  Ed. 

- ,  and  - ,  1973:  Coastal  upwelling  on  a  0-plane.  j  Phys. 

Oceanogr .  3,  16-32. 

Lai,  D.  Y..  1984:  Mean  flow  and  vanabilities  in  the  Deep  Western 
Boundary  Current.  J  Phys  Oceanogr  .  14,  1488-1498. 
Luyten.  J.  R..  J.  Pedlosky  and  H  Stommel.  1983:  The  ventilated 
thermocline.  J  Phys  Oceanogr.  13,  292-309 
Marsh.  J.  G..  R.  E.  Cheney.  J.  J.  McCarthy  and  T.  V,  Martin.  1984: 
Regional  mean  sea  surface  based  upon  GEOS-3  and  Seasat  al¬ 
timeter  data.  Mar  Geod .  8,  385-402. 

McCartney.  M.  S.,  I..  V.  Worthington  and  M.  E  Raymer.  1980 
Anomalous  water  mass  distribution  at  55  W  in  the  North  At¬ 
lantic  in  1977  J  Mar  Res  38,  147-172 
Mellor.  G.  1...  C  R.  Mcchoso  and  E.  Keto.  1982:  A  diagnostic  cal¬ 
culation  of  the  general  circulation  of  the  Atlantic  Ocean  Deep 
Sea  Res.  29,  1 171-1 192 

Ntilcr.  P  P..  and  C  J.  Kohlinsky.  1985:  A  local  time  dependent 
Sverdrup  balance  in  the  Northeast  Pacific.  Science.  229,  754- 
756. 

Olbcrs.  D  J  .  M  Wenzel,  and  J  Willcbrand.  1985:  T  he  inference  of 
North  Atlantic  circulation  from  climatological  hydrographic 
data  Res  (ieoplns  .  23,  313-356 
Orlanski.  I  .  1976  A  simple  boundary  condition  for  unbounded  hv 
perbohe  flows.  J  Comput  Phvs  .  21.  251-269 
Richardson.  P  I  .  19X1  Gulf  Stream  trajectories  measured  with  tree- 
drifting  buoys  ./  Pins  (heanogr  11,990-1010 
.  19X3  Eddy  kinetic  energy  in  the  North  Atlantic  from  surface 
drifters  .1  (ieoplns  Res  88,  4J55-4367 
.  1985.  Average  velocitv  and  transport  ol  the  Gulf  Stream  near 
55  W  J  \Utr  Res  43.  XT- 1  I  I 

Richardson.  W.  S  .  W  J  Schnutz  and  P  P  Nuler.  19n9  I  he  velocity 
structure  of  the  Florida  Current  from  the  Straits  of  Honda  to 
Cape  Fear  Deep  Sea  Res  !6(  Suppl  I,  225-2 3 1 
Riser.  S  (  II  Freeland  and  II  I  Rossby.  I97X:  Mesoscale  motions 
near  the  deep  western  boundary  current  of  the  North  Atlantic 
Deep  Sea  P.  -.  25,  1  1 1191 

Schmitz.  W  J  ,  1 9X0  Weakly  depth-dependent  segments  of  the  North 
Atlantic  circulation  J  Mar  Res  38,  1 1  I  - 1 33 
.  1984:  Abvssal  edds  kinetic  cncrgs  in  the  North  Atlantic.  J 
Mar  Res  42.  509-536 


L 


814 


JOURNAL  OF  PHYSICAL  OCEANOGRAPHY 


Volume  19 


- ,and  W.  R,  Holland,  1982:  A  preliminary  comparison  of  selected 

numerical  eddy-resolving  general  circulation  experiments  with 
observations.  J.  Mar.  Res..  40,  75-1 17. 

- ,  and - ,  1986:  Observed  and  modelled  mesoscalc  variability 

near  the  Gulf  Stream  and  Kuroshio  Extension.  J.  Geophys.  Res., 
91,  9624-9638. 

- ,  W.  R.  Holland  and  J.  F.  Price,  1983:  Mid-latitude  mesoscalc 

variability.  Rev.  Geophys.  Space  Phys.,  21,1 109-1 1 19, 

Stommcl,  H.  M.,  1958:  The  abyssal  circulation,  Deep-Sea  Res.,  5, 
80-82. 

Swallow,  J.  C.,  and  L.  V.  Worthington,  1961:  An  observation  of  a 
deep  counter-current  in  the  western  North  Atlantic.  Deep-Sea 
Res..  8,1-19 . 

Thompson,  J.  D.,  1986:  Altimeter  data  and  geoid  error  in  mesoscalc 
ocean  prediction:  some  results  from  a  primitive  equation  model. 
J.  Geophys.  Res..  91,  2401-2417. 

- ,  and  H.  E.  Hurlburt,  1982:  A  numerical  study  of  the  influence 

of  the  New  England  Seamount  Chain  on  the  Gulf  Stream:  Pre¬ 
liminary  results.  Proc.  Workshop  on  Gti(f  Stream  Structure  and 
Variability.  Univ.  of  N.C.,  Chapel  Hill,  J.  M.  Bane.  Jr..  Ed., 
Office  of  Naval  Research,  346-362. 


Vun  Wyckhouse,  R.  J.,  1979:  Synthetic  Bathymetric  Profiling  System. 
Volume  I:  Data  sources  and  data  preparation.  Naval  Ocean  Re¬ 
search  and  Development  Activity,  NSTI,  Ms.,  NORDA  Tech. 
Note  35, 

Vcronis,  G.,  1981:  Dynamics  of  large-scale  ocean  circulation.  Evo¬ 
lution  of  Physical  Oceanography,  B.  A.  Warren  and  C.  Wumcto. 
Eds.,  The  MIT  Press,  157  pp. 

Wallcraft,  A.  J.,  1980:  Capacity  matrix  techniques.  Ph.D.  Thesis. 
University  of  London,  Imperial  College,  267  pp. 

- ,  and  J.  D.  Thompson,  1984:  Ocean  modelling  and  drifters. 

1984  D'ifliny  Huoy  Workshop  Proceedings.  NaL  Space  Tech. 
Lab ,  Miss.  Mar.  Tech.  Soc.,  81-98. 

Worthington.  L.  V„  1962:  Evidence  for  a  two  gyre  circulation  in  the 
North  Atlantic.  Deep-Sea  Res..  9,  51-67. 

Wunsch,  C.,  and  B.  Grant,  1982:  Towards  the  genera!  circulation  of 
the  North  Atlantic  Ocean.  Progress  in  Oceanography.  Vol.  1 1. 
Pcrgnmon,  1-59. 

Wyrtki,  K„  L.  Magaard  and  J.  llagar,  1976:  Eddy  energy  in  the 
ocean.  ./.  Geophys.  Res.,  81,  2641-2646. 


